Gapless Phase in the XYZ Spin-1/2 Chain in a Magnetic Field and the Quantum 
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We study the XYZ spin-1/2 chain placed in a magnetic field pointing along the a;-axis. We use 
bosonization and a renormalization group analysis to show that the model has a non-trivial fixed 
point at a certain value of the XY anisotropy a and the magnetic field h. Hence, there is a line of 
critical points in the plane (a, h) on which the system is gapless, even though the Hamiltonian has 
no continuous symmetry. The quantum critical line separates a gapped commensurate phase from a 
gapped incommensurate phase. Our study explains why the floating phase of the axial next-nearest 
neighbor Ising (ANNNI) chain in a transverse magnetic field is only a line, as shown by recent 
numerical studies. 
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PACS number: 75.10.Jm, 64.70.Rh 

One-dimensional quantum spin systems have been 
studied extensively ever since the isotropic spin-1/2 chain 
was solved exactly by Bethe. Baxter later used the Bethe 
ansatz to solve the fully anisotropic (XYZ) spin-1/2 
chain in the absence of a magnetic field [|lj; the prob- 
lem has not been analytically solved in the presence of 
a magnetic field. Experimentally, quantum spin chains 
and ladders are known to exhibit a wide range of un- 
usual properties, including both gapless phases with a 
power-law decay of the two-spin correlations and gapped 
phases with an exponential decay There are also 

two-dimensional classical statistical mechanics systems 
(such as the axial next-nearest neighbor Ising (ANNNI) 
model) whose finite temperature properties can be under- 
stood by studying an equivalent quantum spin-1/2 chain 
in a magnetic field. The ANNNI model has been stud- 
ied by several techniques, and it was believed for a long 
time to have a floating phase of finite width in which 
the system is gapless [Q. However, a recent large-scale 
numerical study has shown that the width of this phase 
seems to be zero within numerical errors 1^. 

Amongst the powerful analytical methods now avail- 
able for studying quantum spin-1/2 chains is the tech- 
nique of bosonization Recently, the XX Z chain in 
a magnetic field |^ and the quantum ANNNI model 0] 
have been studied using bosonization. In this work, we 
will study the fully anisotropic XYZ model in a magnetic 
field pointing along the x-axis. For small values of the 
XY anisotropy a and the magnetic field h, we will show 
that there is a non-trivial fixed point (FP) of the renor- 
malization group (RG) in the (a, h) plane; the system is 
gapless on a quantum critical line of points which flow to 
this FP. The floating phase of the ANNNI model wiU be 
shown to be a special case of our results corresponding 
to A = 0. The gapless line is somewhat unusual because 
the XY anisotropy and the magnetic field both break 



the continuous symmetry of rotations in the x — y plane. 
We will provide a physical understanding of this line by 
going to the classical (large S) limit of the model. 

We consider the Hamiltonian defined on a chain of sites 



H = 



[{1 + a) SlSl^^ + (1 - a) SlSl^^ 



hS^. 



(1) 



We will assume that the XY anisotropy a and the zz cou- 
pling A satisfy —1 < a, A < 1. We can assume without 
loss of generality that the magnetic field strength h> 0. 
For a — h — 0, the model is symmetric under rotations 
in the x — y plane and is gapless. The low-energy and 
long-wavelength modes of the system are described by 
the bosonic Hamiltonian IHl 



Ho - 



dx [ {d^ef + {d,(j>f 



(2) 



where v is the velocity of the low-energy excitations 
(which have the dispersion oj — v\k\); z; is a function of A. 
(The continuous space variable x and the site label n are 
related through x — nd, where d is the lattice spacing). 
The bosonic theory contains another parameter called K 
which is related to A by 



K ^ 



7r + 2sin"^(A) 



(3) 



K takes the values 1 and 1/2 for A = (which describes 
noninteracting spinless fermions) and 1 (the isotropic an- 
tiferromagnet) respectively; as A ^ —1, K ^ oo. We 
thus have 1/2 < < oo. 

The spin-density operators can be written as |p| 



^ \17 + (-l)"ciC0s(2V7ri^,/)) , 
S^i = [c2Cos{2V^(b) + (-1)"C3] cos(,&) , (4) 

V A 



1 



where the c; are constants given in Ref. 10. The XY 
anisotropy term is given by 



where C4 is another constant. 

For convenience, let us define the three operators 



01 = cos{2VTTK(t))cosU —e) , 

V K 

02 = cos{2, ^9) , and O3 = cos{aV^4>) . (6) 

V K 

Their scaling dimensions are given by K + 1/AK, 1/K 
and AK respectively. Using Eqs. the terms corre- 

sponding to a and h in Eq. dl]) can be written as 



Ha + Hh — 



dx [ ac402 — hc20i 



(7) 



where we have dropped rapidly varying terms propor- 
tional to (—1)" since they will average to zero in the 
continuum limit. (We will henceforth absorb the factors 
C4 (C2) in the definitions of a (h)). We will now study 
how the parameters a and h flow under RG. 

The operators in Eqs. @ are related to each other 
through the operator product expansion; the RG equa- 
tions for their cocfhcients will therefore be coupled to 
each other [|ll| . In our model, this can be derived as fol- 
lows. Given two operators Ai — exp(iai(/) -I- iPi9) and 
A2 — exp(iQ!20 -l- we write the fields (j) and 9 as 

the sum of slow fields (with wave numbers |fc| < Ae^'") 
and fast fields (with wave numbers Ae~^' < |fc| < A), 
where A is the momentum cut-off of the theory, and dl is 
the change in the logarithm of the length scale. Integrat- 
ing out the fast fields shows that the product of Ai and 
A2 at the same space-time point gives a third operator 
A3 = e*("i+"2)0-Hi(/3i-H/32)e ^i^j^ prefactor which can be 

schematically written as 



A1A2 



-(Qia2+/3i/32)<i//27r a 



(8) 



If Xi{l) denote the coeffients of the operators Ai in an 
effective Hamiltonian, then the RG expression for dX^ / dl 
will contain the term (aia2 -I- /3i/32)AiA2/27r. Using this, 
we find that if the three operators in (^ have coefficients 
h, a and b respectively, then the RG equations are 



da 

'dl 

db 

di 

dl 4 



(2 - —)a - (2K - —)h'^ 
^ k' ^ 2K' 

(2-4K)6 + {2K-^)h^ 



(9) 



where we have absorbed some factors involving v in the 
variables a, b and h. (Wc will ignore the RG equation 



for V here). It will turn out that K renormalizes very 
little in the regime of RG flows that we will be concerned 
with. Eqs. ^ have appeared earlier in the context of 
some other problems [||J^]. However, the last two terms 
in the expression for dh/dl were not presented in Ref. 
8; these two terms turn out to be crucial for what fol- 
lows. Note that Eqs. (^ are invariant under the duality 
transformation if <-> 1/ AK and a *^ b. 

Let us now consider the fixed points of Eqs. (||). For 
any value of i^T = K* , a trivial FP is {a*,b*,h*) — 
(0,0,0). Remarkably, it turns out that there is a non- 
trivial FP for any value of K* lying in the range 1/2 < 
K* < 1 + V3/2 ; we will henceforth restrict our attention 
to this range of values. (The upper bound on K* comes 
from the condition 2-i\r* -l/AK* > 0). The non-trivial 
FP is given by 



h* = 



.y2K*{2-K* - 1/4K*) 
{2K* + 1) 



a* = {K* + -)h*^ , and b* 



2K* 



(10) 



The system is gapless at this FP as well as at all points 
which fiow to this FP. One might object that Eqs. (||) 
can only be trusted if a, b and h are not too large, other- 
wise one should go to higher orders. We note that the FP 
approaches the origin as K* ^ 1 + V3/2 ~ 1.866; this 
corresponds to the zz coupling A = — sin[7r(\/3 — 3/2)] ~ 
—0.666. Thus the RG equations can certainly be trusted 
for K* close to 1.866. For K* = 1, the FP is at 
{a*,b*,h*) = (1/4,1/8,1/V^). 

We have numerically studied the RG flows given by 
Eqs. (^ for various starting values of {K,a,b,h). Since 
the Hamiltonian in (|^) does not contain the operator O3, 
we set & = initially. We take a and h to be very small 
initially, and see which set of values flows to a non-trivial 
FP. For instance, starting with K — 1, b — 0, and a, h 
very small, we find that there is a line of points which fiow 
to aFP at {K\a*,b\h*) = (1.020,0.246,0.122,0.404). 
This line in the (a, h) plane is shown in Fig. 1. We see 
that K changes very little during this flow; if we start 
with a larger value of K initially, then it changes even 
less as we go to the non-trivial FP. It is therefore not a 
bad approximation to ignore the fiow of K completely. 

We can characterize the set of points (a, h) lying close 
to the origin which flow to the non-trivial FP. Numeri- 
cally, we find that there is a unique fiow line in the (a, h) 
plane for each starting value of K and 6 = 0, provided 
that a, h are very small initially. This means that a(l) and 
h{l) given by Eqs. (H) must follow the same line regard- 
less of the starting values of a,h. From Eqs. (||), we see 
that if h « then h{l) ~ h{0) exp(2 - K - l/AK)l 
while a{l) ~ exp(2 — 1/K)l. Hence h must initially scale 
with a as 



(11) 



as we have numerically verified for K — 1. However, Eq. 
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FIG. 1. RG flow diagram in the (a, h) plane. The solid line 
shows the set of points which flow to the FP at a* — 0.246, 
h* — 0.404 marked by an asterisk. The dotted lines show the 
RG flows in the gapped phases A and B (see text). 



(y) is only true if (2 - - l/4:K)/{2 -l/K) > 1/2, 
i.e., if K < (1 + V2)/2 = 1.207. For K > 1.207 (i.e., 
A < —0.266), the initial scaling form is given by /i ~ a^^^. 

We now examine the stability of small perturbations 
away from the fixed points. The trivial FP at the 
origin has two unstable directions (a and h), one sta- 
ble direction (b) and one marginal direction (K). The 
non-trivial FP has two stable directions, one unstable 
direction and a marginal direction (which corresponds 
to changing K* and simultaneously a*, b* and h* to 
maintain the relations in Eqs. (p^)). The presence 
of two stable directions implies that there is a two- 
dimensional surface of points (in the space of param- 
eters (a, b, h)) which flows to this FP; the system is 
gapless on that surface. A perturbation in the unsta- 
ble direction produces a gap in the spectrum. For in- 
stance, at the FP {K\a*,b\h*) = (1,1/4,1/8,1/^6), 
the four RG eigenvalues are given by 1.273 (unstable), 
(marginal), and -1.152 ± 1.067? (both stable). The 
positive eigenvalue corresponds to an unstable direction 
given by {SK,6a,db,6h) = (5a(0.113, 1, -0.092, -0.239). 
A small perturbation of size 6a in that direction will 
produce a gap in the spectrum which scales as AE 
|ju|i/i.273 _ |ju|0.786. -(-j-^g correlation length is then given 

Fig. 1 shows that the set of points which do not flow to 
the non-trivial FP belong to either region A or region B. 
These regions can be reached from the non-trivial FP by 
moving in the unstable direction, with Sa > for region 
A, and 6a < for region B. In region A, the points flow 
to a = oo; this corresponds to a gapped phase which is 



commensurate since the magnetic field h asymptotically 
flows to zero. In region B, both a and h flow to — oo; this 
is a gapped incommensurate phase. The line separat- 
ing the two regions is a commensurate-incommensurate 
transition line. An order parameter which distinguishes 
between regions A and B is the staggered magnetization 
in the y direction, defined as 



= [ hm (-1)' 



< SqS^ > 



il/2 



(12) 



This is zero in the commensurate phase A. In the in- 
commensurate phase B, it is non-zero, and its scaling 
with the perturbation Sa can be found as follows Q . At 
a = h = 0, the leading term in the long-distance equal- 
time correlation function of is given by 



<sy 



(-1)" 

|„ll/2/f 



(13) 



Hence the scaling dimension of S^^ is 1/4K. In a gapped 
incommensurate phase in which the correlation length 
is much larger than the lattice spacing, my will therefore 
scale with the gap as niy ~ (AE)^^'^-^ . If we assume that 
the scaling dimension of S*-^ at the non-trivial FP remains 
close to l/4:K, then the numerical result quoted in the 
previous paragraph for K — 1 implies that my ~ IJaj*^'^^^ 
for 6a small and negative. 

We will now apply our results to the one-dimensional 
spin- 1/2 quantum ANNNI model with nearest 

neighbor ferromagnetic and next-nearest neighbor anti- 
ferromagnetic Ising interactions and a transverse mag- 
netic field. The Hamiltonian is given by 



E 



r 



[ - 2J,S:S^^, + J2StSt+, + -5f ] , (14) 



where Ji, J2 > 0; we can assume without loss of gen- 
erality that r > 0. The quantum Hamiltonian in (|l^ ) 
is related to the transfer matrix of the two-dimensional 
classical ANNNI model; the finite temperature critical 
points of the latter are related to the ground state quan- 
tum critical points of (p^), with the temperature T being 
related to the magnetic field F. Earlier studies showed 
that the model has a fioating phase of finite width which 
is gapless and incommensurate A recent bosoniza- 
tion study reached the same conclusions 0. However, 
a more recent numerical study involving very large two- 
dimensional classical systems at finite temperature con- 
cluded that the floating phase seems to have zero width 
We can now understand this as follows. Consider a 
Hamiltonian which is dual to ( p^ for spin- 1/2; this will 
turn out to be a special case of our earlier model. The 
dual Hamiltonian is given by P,p2[ 



Hd = J^TnT^l+l + ^nT'n+l - JiT: ] , (15) 

n 

where T" are the spin-1/2 operators dual to S". Let us 
scale this Hamiltonian by a factor to make J2 + F = 2. 
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Then this has the same form as in Eq. with J2 = 
1 + a, T = 1 — a, Ji = h, and A = 0. Hence it follows that 
the quantum ANNNI model has a line of points in the 
( J2 — r, Ji) plane on which the system is gapless. (From 
Eq. (pi]), we see that the shape of this line is given by 
Ji ~ (J2 - r)3/4 as Ji ^ 0). This line can be identified 
with the floating phase. Our results for the width of the 
floating agree completely with those obtained for the two- 
dimensional classical ANNNI model in Ref. 5. We should 
point out, however, that Ref. 5 shows the phase transi- 
tion to be of the Kosterlitz-Thouless type (with ^ diverg- 
ing exponentially) from the high-temperature side (i.e., 
from region B in Fig. 1), and is of the commensurate- 
incommensurate type (with ^ diverging as a power-law) 
from the low- temperature side (region A in Fig. 1). Our 
analysis indicates that it should be of the commensurate- 
incommensurate type from both sides. This discrepancy 
remains to be resolved. One comment to be made here is 
that Ref. 5 has studied two values of J2/J1 (0.6 and 0.8) 
which are not large; our RG results may be expected to 
be valid if a, h are small, i.e., if J2/ Ji is large. 

Finally, we would like to provide a physical picture 
of the gapless line in the (a, h) plane by looking at the 
classical limit of Eq. (|l]) . Consider the Hamiltonian 

n 

+ As:,s:,+, - hss: ] , (le) 

where the spins satisfy = S{S + 1), and we are inter- 
ested in the classical limit S" — > 00 |12| ]. Let us assume 
that the zz coupling is smaller in magnitude than the xx 
and yy couplings. Then the classical ground state of (|l^ ) 
is given by a configuration in which all the spins lie in the 
X — y plane, with the spins on odd and even numbered 
sites pointing respectively at an angle of ai and —a2 to 
the X-axis. The ground state energy per site is 

e(ai,Q;2) — S'^ [— — (cos ai -|- cos 02) 

-I- cos(q;i -I- q;2) + acos(Q;i — 02)]. (17) 

Minimizing this with respect to ai and a2, we discover 
that there is a special line given by = 16a on which 
all solutions of the equation 



hcos{ 



4a cos( 



ai — a2 - 



(18) 



give the same ground state energy per site, cq ^(1 + 
a)S^. The solutions of Eq. ( |l8| ) range from ai — a2 — 
cos~^(4a//i) to ai — TT,a2 = (or vice versa); in the 
ground state phase diagram of the ANNNI model, these 
two configurations correspond to the paramagnetic in- 
commensurate phase and the antiphase respectively 
Classically, we see that the line = 16a separates these 
two phases; we therefore identify this line with the float- 
ing phase The fact that there is a one-parameter 



set of ground states (characterized by, say, the value of 
ai in the solution of (^8|)) which are all degenerate for 

= 16a means that there is a gapless mode in the exci- 
tations; we can derive this explicitly by going to the next 
order in a 1/S expansion. Hence the model is gapless 
on that line. This provides some understanding of why 
one may expect such a gapless line in the spin-1/2 model 
also. Note however that the bosonization analysis gives 
the scaling form in (^Tj) for h versus a; this agrees with 
the classical form only if A < —0.266. 

To summarize, we have shown that the XYZ spin-1/2 
chain in a magnetic field exhibits a gapless phase on a 
particular line. It may be interesting to use numerical 
techniques like the density-matrix renormalization group 
method |l^] to study various ground state properties of 
this model, in particular the behavior of order parameters 
such as the one defined in Eq. (|l^). The RG equations 
studied in this paper appear in other strongly correlated 
systems, such as the problem of two spinless Tomonaga- 
Luttinger chains with both one- and two-particle inter- 
chain hoppings H] , and one-dimensional conductors with 
spin- anisotropic electron interactions The gapless 
phase may therefore also appear in other systems. 
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